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I. INTRODUCTION 

Optical conductivity data on <j{lo) vs cj contains important information on electron dynamics. In general however the 
relationship between a{ui) and the electron self-energy is rather complicated. Under such circumstances the reduction, 
if possible, of the complete expression for conductivity, to somewhat more approximate but simpler analytic form can 
be valuable. It can help our understanding of the basic physics as well as provide experimentalists with a simpler 
basis for the analysis of data. 

In the specific case of an electron-phonon system, the analytic formula provided by Allen [l| has proved to be 
valuable. It relates the measured optical scattering rate through a simple integral to the underlying electron-phonon 
spectral function a 2 F(uj) which is, in the end, the fundamental quantity of interest. A generalization of Allen's 
formula to finite temperature was provided by Shulga et. al. Q ■ In this paper we want to extend this previous work 
to the case when the underlying electronic density of states is energy dependent rather than constant. A motivation 
for this extension is to provide guidance in the interpretation and analysis of optical data in the pseudogap regime of 
the cuprates. In fact the formula derived herein has already been used in the experimental work of Hwang et. al. Q 
on ortho-II YBCO. 

II. PRELIMINARIES 



The Drude formula for optical conductivity a(u) = ai(uj) + io2(w) 

a(w,T) = ^ — - (1) 

can be extended (see 0, Q and Refs. therein) to include a frequency dependent scattering rate 

UJp 1 

T) = 4^ 1/t(lu, T) - ko[l + fj] ' (2) 

where 1/t(uj,T) is the frequency-dependent optical scattering rate and X(uj,T) is the optical mass enhancement. For 
a spherical Fermi surface the plasma frequency uj 2 — Anne 2 /m, where n is the free-carrier density and m is the carrier 
mass. 

One can solve Eq. Q for 1/t(oj) and 1 + \(oj) in terms of the optical conductivity found from experiment, 
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and 



u 2 1 

1 + A(w) = -T^-Im 
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The plasma frequency can be extracted from the experimental data using the sum rule a\{ui)djuj = uip/8. Although 
the representations of experimental data using o~i(u>), 0*2(0; ) and 1/t(u>) with 1 + A (a;) are formally equivalent, it has 
become rather popular to discuss the pseudogap behavior in high-temperature superconductors (HTSC) using the 
language of the optical scattering rate and mass enhancement. For example, a drop of l/r(w) extracted from in-plane 
optical conductivity measurements in HTSC which is observed below a certain frequency for the temperatures T < T* b 
is associated with the above-mentioned pseudogap 0. 

Another advantage of 1/t(uj) is that in an electron-phonon system it is related to the electron-phonon interaction 
spectral density, a 2 F(uj). For example, there is an approximate relationship 6] 



t(uj) 



(5) 



which is valid at T = in the normal state. Note that one can consider a general electron-boson interaction function, 
because instead of a 2 F(uj) one can, for instance, introduce electron-spin excitation spectral density, I 2 x{^) [3) Thus 
below we imply this more general case, but preserve the historical notation a 2 F(uj). 

Despite its "magic" simplicity, formula JSJ) works rather well (see e.g. 0, but unfortunately its practical 
applications are limited because the experimental data must be (ambiguously) smoothed "by hand" before the second 
derivative can be taken. This problem was very recently circumvented in Ref. by solving the corresponding 
integral equations for a 2 F(uj). For example, instead of differentiating uj/t(uj), one can extract a 2 F(uj) using the 
well-known result of Allen [l| 
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The last expression was derived using the second order perturbation theory and in fact, Eqs. JSJ and (JSj are equivalent 
under the condition 1/t(uj = 0) = 0. 

A finite temperature generalization of Eq. © was derived in Ref. [2| using the Kubo formula 
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In fact, there is no finite temperature equivalent of a"magic" equation iJSJ for Eq. (J7J, so that for finite T's a 2 F can 
only be found by inversion of the last equation. The numerical method of inversion of Eq. Q, its limitations and the 
resulting doping and temperature dependences of the bosonic spectral function, a 2 F{uj) in several families of HTSC 
were investigated in detail in Ref. [l0| . 

There is, however, an important assumption used in deriving both the T = equation © and its finite temperature 
extension l(7)l. viz. the electronic density of states (DOS), N{e) is taken as a constant in the vicinity of the Fermi 
level. As shown in Ref. [ljj when Eq. Q is used for analyzing the experimental data in the pseudogap state of 
HTSC, the resulting a 2 F contains nonphysical negative values. This problem originates from the fact that the 
above-mentioned assumption N(e) — const is definitely strongly violated in the pseudogap state. The essence of the 
pseudogap phenomenon is that N(e) becomes a nontrivial function of energy and the form of N(e) depends strongly 
on the temperature and doping. 

Therefore it would be very useful to have a generalization of Eqs. I© and 0) valid for iV(e) 7^ const. Interestingly, 
such a generalization of the zero temperature expression JBJ was already done twenty years ago by Mitrovic and 
Fiorucci in relation to A15 compounds 
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Eq. ((SJ) was derived using the method of Ref. [l| and it is easy to see that for N(e) — const it reduces to Eq. ©. 

The purpose of the present work is to obtain a generalization of the finite temperature expression J7J| valid for 
N(e) ^ const. In contrast to Refs. we base our derivation on the Kubo formula which turns out to be more 

useful for considering the T ^ case. 

We begin by presenting in Sec. IIIII the expression for optical conductivity c(uj) in terms of frequency dependent 
self-energy S(w) for the case N(e) ^ const. The corresponding expression for T,(uj) is obtained in Sec. II VI for the 
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case of non-constant quasiparticle DOS N(e). The difference between the usual DOS N(e) and the quasiparticle DOS 
N(e) is pointed out. In Sec. [V] we present the relationship between the optical scattering rate, l/r(w), electron-boson 
interaction function a 2 F{uj) and the quasiparticle DOS N(e). The frequency dependent impurity scattering rate 
l/ri mp (w) is considered in Sec. I VII In the Discussion, Sec. IVIll we illustrate that the pseudogap opening results in 
the decrease of 1/t(u>) and consider possible applications of our results. 



III. OPTICAL CONDUCTIVITY FOR N(£) const 
We begin with the Kubo formula for electrical conductivity ^3] 

°a M = — [nj - (w + io)] , (9) 

where r,j is the diamagnetic term and Ily (u) is the retarded correlation function obtained by analytical continuation 
[n^(a;) = Hij(iQ m — > uu + iO)] of the imaginary time expression 

n y (^ m ) = i [^ dTe m ™ T (ji(q_ = 0,^(0,0)), O m = ^. (10) 
•'Jo P 

Here ji(q, t) is the Fourier transform of the paramagnetic electric current density operator, V is the volume of the 
system, T = 1/(3 and for the case of parabolic band Ty = uip/4irSij — e 2 n8 i j/m. The electrical conductivity @ 
consists of a regular part 

o^(a;) = -i[n«H-ng(o)] (11) 

and singular part, related to the superconducting condensate, 

4 c H = -^N- n *(°)]- (12) 



Since in what follows we restrict ourselves to considering the normal state of the isotropic two dimensional system, 
we suppress the subscript reg and consider o{uj) = a™Jf(w) = er™ s (ti>) deriving, for example, H xx (u>). 
Neglecting vertex corrections, the calculation of <j(uj) reduces to evaluation of the bubble diagram 

°° f d 2 k 

Hij(i£l m ) = -2e 2 T ^2 / j^-T2 v FiQ*)vFjQi)G(iu) n +iO m ,k)G(iw„,k), (13) 

where i>Fi(k) = <9£(k)/<9fci is the Fermi velocity and 

G(iwn, k) = , 1 (14) 

iLu n -f(k) - £(iw n ,k) 

is the fermionic Green's function with the self-energy S(iu; n ,k) and u„ = 7r(2n + l)//3. Using the spectral represen- 
tation for the Green's function l|14(l . one can easily sum over fermionic Matsubara frequencies in Eq. (| 1 31) 

n xa (ffi m ) = ~2e 2 ^ d£ [ ^L S (Z - £(k))v% x (k) r dio 1 r du 2 UF ^ ] - A(m , k) A(o; 2 , k) , (15) 



where the spectral function 



A(w,k) = --lmG R (u + iO,k), (16) 

7T 



hif(w) = l/(eP u + 1) is the Fermi distribution, and we inserted the integral over £ which is equal to 1. Since we 
are interested in the complex conductivity <j(lo) we cannot remove one of the integrations over w in i|15[l by taking 
ImII(fjj + iO), but one of the integrations can be done by again using the spectral representation for retarded (advanced) 
Green's function 

G r ,a{v,k)= du 2 t - n • ( 17 
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Then we obtain 



/oo />oo 
du'n F (u') d4N(0v 2 Fx (0A(u',0[GR^' + ^,0 + GA(u>'-u > ,0}, (18) 
-oo J — oo 



where to isolate the effects of the energy dependence of the single-spin-band DOS, 



n(o = J -cm as) 

the velocity v Fx (t;) is defined as (see Refs. [lllll3$ l 

v Fx (0 = J -0^v Fx (k)S(( - e(k)). (20) 

Writing Eq. (|18|l we also assumed that A{ui,k) and, accordingly, Gij^(a;,k) to be dependent only on £(k). 

Now we must make two important assumptions: The first one is quite common and states that the self-energy 
£(za;,£) does not depend on £, so that the whole dependence of Gb.,a(^>,S,) is contained in the free-electron dispersion 
£(k). The second assumption is that the energy dependence of the square of the plasma frequency 



! -2e'N(t)v Fx (0 (21) 
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can be ignored as compared to the dependence of N(£) in the vicinity of £ = 0, so that in Eq. (|18H we can replace 
^piO by w p(£ = 0)- The validity of this approximation for A15 compounds was discussed in Ref. [llj and here we 
will assume that it is also valid for HTSC. In this respect, the first assumption can be considered as a statement that 
£(o;, £) is approximated by E(w, £ = 0) [lq. After these two assumptions are made we can integrate over £ and finally 
arrive at the following representation for the optical conductivity (see e.g. Refs. pi ITil HHUi! ) 

cr(w) = -^-/ de[n F (e) - n F (e + u)} — — — -, (22) 

iiruj^ l w v n w + z/T imp + S*(e) - E(e + cj)' v ; 

where E(e) is the retarded self-energy on the real axes and S*(e) its complex conjugate. In Eq. we also included 
the electron-impurity scattering rate, l/ri mp , its frequency dependence will be considered in Sec. I VII One can check 
that for T — Eq. (|2*2)) reduces to the expression for cr(uo) written in Refs. [EIEI3- 

To investigate the effect of electron-boson interaction on l/r(w) we must express Re[l/<7(w)] in terms of the self- 
energy £ in the simplest possible form. It can be anticipated if one substitutes Eq. (1221 in Eq. J2J and rewrites it as 
follows: 

1 f 1 f°° 1 I -1 
-— =wW — — / de n F (e) -n F {e + u)]- — — — r — — — — \ (23) 

Now expanding the denominator of Eq. (|23|) . doing the integration and then "de-expanding" the result (see Refs. 
we obtain the following approximate representation 



1 1 1 



de[n F {e) - n F (e + u>)] Im[S(e + uj) - S*(e)]. (24) 



t(u>) T imp u 

In deriving Eq. (|24|l we used the assumption |E(e + u>) — S*(e)| <C \oj + «/Ti mp | to expand and then assumed that 



de [n F {e) — n F (e + lu)\ |S(e + ui) — S*(e)| <C u)\oj + i/ r im P | 



to "de-expand" . Based on these inequalities one would expect that all results that follow from Eq. I|24|l are valid 
only for large oj. Nevertheless, the numerical comparison of the results obtained by the direct substitution of the 
self-energy (|32H in Eq. H22fl with 1/t(cj) computed using Eq. Q in Ref. |2j shows that the last equation is valid for 
a much wider range of the frequency u>. This is in spite of the fact that the derivation of Eq. J7|) is based on the 
approximate Eq. I|24|l . Finally we remind that since we did not include vertex corrections, 1/t(oj) is expressed in terms 
of the usual self-energies instead of transport "self-energies" discussed in Ref. Accordingly in Sec. [3 the optical 
scattering rated l/r(w) will be expressed in terms of the tunneling a 2 F instead of the transport a 2 r F considered by 
Allen in Ref. Q. 
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IV. SELF-ENERGY FOR N(£) / const 

Now we consider the influence of a nonconstant DOS on the usual relationship between the self-energy T,(u>) and 
the electron-boson interaction function a 2 F. We begin with the well-known expression (see e.g. Refs. [l7l fl8| ) 

E(iw n )=T / d ^jm dCla 2 F(Cl) Q2 + Q2 G(iw n - iti m ,Z) (25) 

rn — — 00 ~ 00 V / m 

written for the case when a 2 F(Q) does not depend on the electron energy. In Eq. (|25f) 

G(iu n , = j^J 7S2 § ti - e(k))G(t Wn , k), (26) 

the DOS iV(£) is defined by Eq. I(19() and its energy dependence usually is also neglected. Our goal is, however, to 
retain N(£) and consider the influence of iV(£) 7^ const on E(w) and, accordingly, on a(u>). Again using the spectral 
representation for the Green's function G(iu; n ,k), one can easily sum over Matsubara frequencies in Eq. ((25(1 and 
obtain 

E(iw„) = J ^SSy dw'J^ dCla 2 F(Cl) ^~hnG R (u + iO,Z)^I(iu n ,n,L>'), (27) 



where llj 



, ,v _ n B (n) + l-n F (a; / ) n B (») +n F (^) 
'K. S! . W )- ; p; ; 1 : — pi — l^oj 



with the Bose distribution ns(0) = l/(e' 3 ° — 1). 
Let us now consider the quantity 

N(u) = -- f dZN(£)TmG R (u + iO,£) (29) 

77 J-oo 

which enters Eq. I|27|l . Using Eq. 1(26(1 and the definition of the spectral function one can easily check that 

N(w) = J -0^A(u,k), (30) 

viz. this quantity represents the fully dressed quasiparticle DOS which could contain a pseudogap that has its origin 
in correlation effects. Finally making an analytical continuation iw n — > lj + iO and taking the imaginary part of E(w) 
we obtain 



Im£H = -tt jH dCla 2 F(Cl) j ^_"l[„ B (n) + 1 - n F (uj - CI)] + ^±^[n s (0) + n F {w + CI)] 



1 N(0) v /J JV(0) 

It is easy to see that for N(u>) — N(0) — const the previous equation reduces to a more familiar expression 0,^1 



(31) 



p 00 

ImEfw) = -- dCla 2 F{Cl) 
2 Jo 



CI LO + CI U! — CI 

2 coth tanh h tanh 



2T 2T 2T 



(32) 



OPTICAL SCATTERING RATE: BOSON CONTRIBUTION 



Substituting the self-energy 131(1 in the expression 124(1 and doing simple replacements of the variables, we arrive 
at the main result of the present paper, 



t(u) u J 



dCla 2 F(Cl) / de 



N(e-Cl) N(-e + Cl) 



N(0) 



N(0) 



[n B (Cl) + 71f(CI - e)][n F (e — u) — n F (e + u)], (33) 



which establishes a link between the optical scattering rate 1/t(oj), the electron-boson interaction function a 2 F(Cl) 
and the quasiparticle DOS N(oj). It is important to stress that in contrast to N(u>), the quasiparticle DOS N(ui) 
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which enters Eqs. (|31|l and l|33(l is directly related to the spectral function A(cu, k) measured by ARPES experiments 
[T^ | . Note that in Eq. (|33|) material parameters enter only as the electron-boson spectral density and the fully dressed 
electron DOS. Different mechanisms leading to the same N(oj) are differentiated in the optical scattering rate only 
through the size and shape of a 2 F(fl). 

To understand better the rather lengthy Eq. (|33|l we consider limiting cases where one can establish a link with 
already known results. 

(i) For T = the Bose distribution in Eq. I)33|) drops out and it reduces to Eq. © with the band DOS N(lo) replaced 
by the quasiparticle DOS N(ui). Note that it was pointed out in Ref. that N(uA in Eq. © should be interpreted 
as the quasiparticle DOS, but the golden rule approach of Allen Q used in Ref. ^Jj does not allow to establish this 
fact directly. Comparing Eqs. JSJ and (|33|) one can see that finite temperature brings an essential element, the Bose 
distribution ng(fl). 

(ii) When N(u>) = const it becomes possible to integrate over e in l|33|) . Indeed using the integral 

dzriF(z + a)np(—z — b) = (a — &)rts(a — b) (34) 

we arrive at Eq. (7J obtained by Shulga et al. in Ref. 2]. Obviously, Eq. J7J also follows directly from Eqs. IMl'l) and 

PHi. 

(iii) At temperatures much higher than the boson spectrum upper-energy cutoff, T 3> £l c , expression (|33|l reduces to 



hm * ^r^^r* 



T/n c ^oo t(lo) lu J Q 



N(e-n) N(-e + Q) 
N(0) + N(0) 



[np(e — to) — nf(f + to)]. (35) 



When N(e) = const the last equation can be further simplified || 



hm 4t =4ttT I™ dn^^-. (36) 
t/o c ^oo r(0) J n y ' 

In the case when the electron- phonon interaction is the origin of a 2 F(VL), Eq. I|3f>(l reflects the familiar result that 
the high-temperature electron-phonon contribution to a dc resistivity is linear in temperature. This no longer strictly 
holds for Eq. I|35|l where there is an additional T dependence in the integral over e. 



VI. OPTICAL SCATTERING RATE: CONTRIBUTION OF IMPURITIES 



While Eq. (|33|) represents the electron-boson contribution to the optical scattering rate, the total scattering rate 

1 1 1 



TtotM t(lj) 



(37) 



consists of two parts, viz. the above-mentioned bosonic l/r(w) and that caused by impurities l/ri mp (w). 
There is a simple way of obtaining l/Ti mp (o;) by using the following expression for 

AT imp 7T J 



(38) 



and assuming that the integration of 8(fl) over f2 from to oo gives 1/2. l/ri mp in Eq. I|38|) is the normal state 
impurity scattering rate and it is implied that the limit Q — > in a 2 Fi mp (Q) must be taken before doing the limit 
T — > 0. Substituting Eq. l(55|l in Eq. (|3*3^l we arrive at the expression 



lui 



1 1 

"imo 4w 



de[np(e — lu) — np(e + i 



N(e) N(-e) 
N(0) N(0) 



(39) 



Energy dependence in N(e) implies an energy and temperature dependence in the impurity scattering rate. For T = 
Eq. (SHI reduces to the result of Ref. [ill 



' imp 



l l 

to 



i imp 



de- 



N(e) N(-e) 
AT(0) N(0) 



(40) 
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Note that due to the above-mentioned noncommutativity of the limits O — ^ and T — > in Eq. I|38|) . the last 
expression cannot be obtained simply by substituting Eq. in Eq. JSJl. It is very easy to see that for N(e) = const, 
Eq. (|4*n|) as well as Eq. l|3"§|) reduce to the trivial result Ti mp (aj) = r imp . For a more extensive discussion of the impurity 
scattering problem when the DOS is energy dependent the reader is referred to Ref. [2(| ■ In the fits to data made in 
Ref. 0] residual scattering is small and not important. In that case l/ri mp is estimated to be a few meV while the 
inelastic scattering rate rises to a value larger than 600 meV in the frequency range of interest in the fit. Under such 
circumstances any modulation of 1/Tj mp (ci;) brought about by energy dependence in N(e) is of no consequence when 
fitting to the inelastic scattering part. This would not be so if impurity scattering became large and of the order of 
the inelastic scattering in the important frequency range. 



VII. DISCUSSION 



To illustrate the effect of the opening of the pseudogap on 1 /t(lo) in Fig.^we plot 1 /t{lo) computed with and without 
pseudogap, but do not consider the contribution from impurities. The case without pseudogap which corresponds to 
N(e) = N(0) = const was considered using Eq. 10, while to model a pseudogap we took Q 



N(e) 
N(Q) 



N(0) N(0)\ e 2 

N(0) { N(Q) ) A 2 



0(A-|e|)+0(|e|-A) (41) 



and used Eq. I|3H[I. It is easy to see that the main effect of the opening of the pseudogap [e.g., N(0) < N(0)] is 
to reduce 1/t(cj). This result implies that when there is a pseudogap one cannot use Eq. to estimate a 2 F(fl), 
because the strength of the peaks of a 2 F{Q) would be underestimated. More importantly, the position of these peaks 
would be shifted depending on our assumptions about the absence or presence of the pseudogap. 

Indeed by applying Eq. (|33|) in Ref. k| to the analysis of the experimental data it was demonstrated that there 
is agreement between the sharp bosonic resonance observed in the infrared scattering rate and the properties of the 
(ir, 7r) spin-flip neutron mode. Using the above-mentioned property that the decrease of N(e) due to the opening 
of the pseudogap effectively reduces 1/t(lu), one can choose the position of the bosonic resonance in a 2 F(fl) at 
f2 = 248cm" 1 and thus makes it agree with the frequency of the 31 meV neutron mode, measured by Stock et al. 
[2l|. The use of Eq. (|33[) is crucial for this agreement and the conventional relationship valid for T — and 
iV(e) = const would produce a peak in a 2 F(£l) at £1 — 350cm" 1 . The possibility of achieving a reconciliation between 
the experimental results obtained from neutron scattering and optical conductivity is quite important in developing 
a coherent theoretical description of cuprate superconductors and shows the usefulness of Eq. (|33|l . 

We note that since this kind of analysis involves modeling the form of the pseudogap and its temperature dependence, 
the final results for a 2 F(£l) are definitely not unique. However, one does not expect any qualitative change to the 
model for a 2 F(£l) obtained in Ref. Q which includes a very wide temperature independent background and a sharp 
temperature-dependent peak. As a given sample is studied using several experimental techniques and the quality of 
the data improves, the fit should become more unique. Also it would be interesting to repeat the analysis made in 
Ref. [T^j, but now based on Eq. 13311 derived in this paper, rather than on Eq. J7J) which is valid only for a constant 
density of states. 

Finally we mention that in Ref. |10( data were analyzed not only in the normal or pseudogap state at finite T, 
but also in the superconducting state. This latter analysis is based on a relationship between 1/t(uj) and a 2 F(il) 
derived in Ref. yj for s-wave superconductor at T = Q. In relation to the results obtained in the present paper 
it is worthwhile to ask the following question. Is it possible to distinguish the decrease of \/t(uj) caused by the 
pseudogap and by the superconducting s- or even <i-wave gap? To address this question in detail there is a need to 
generalize the corresponding expression from Ref. Q to d-wave symmetry of the superconducting gap and T =/= 0. 
Such considerations go beyond the scope of the present work. 
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FIG. 1: (Color online) The dependence 1/t(u>) (in arbitrary units) obtained from Eqs. J7J for N(e) — const and 13311 with N(e) 
given by Eq. igTJ. We take the Einstein model for a 2 F(9) = a 2 5(Q, - Q E ) with Q, E = 2A and T = 2.5A. 



